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Beyond the second-order Born approximation, we develop an improved master equation approach 
to quantum transport by virtue of a self- consistent Born approximation. The basic idea is replacing 
the free Green's function in the tunneling self-energy by an effective reduced propagator under the 
Born approximation. We found that the effect of this improvement is remarkable. For instance, 
completely beyond the scope of the conventional second-order master equation, the new approach 
can not only recover the exact result of noninteracting transport under arbitrary voltages, but also 
predict the challenging nonequilibrium Kondo effect. In addition to having an elegant structure, the 
application convenience and accuracy of the proposed scheme, as demonstrated by the examples in 
this work, suggest it a useful tool for quantum transports. 
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I. INTRODUCTION 

Landauer-Biittiker theory and the non-equiHbrium 
Green's function (nGF) approach are widely accepted as 
two typical formalisms for mesoscopic transports ^. 
Besides these two "orthodox" methods, we notice that 
the rate equation approach, an alternative choice, is par- 
ticularly suited for transport though nanostructures with 
few discrete states. Straightforwardly, "classical" rate 
equation can be constructed by phenomenological con- 
siderations Later efforts include its derivation and 
quantum generalization in the context of the resonant 
tunneling system, based on the nGF quantum kinetic 
theory Q, as well as its modification to describe the 
quantum coherence broadly inhered in mesoscopic sys- 
tems ||, 1^. In particular, a microscopic derivation start- 
ing with the many-particle Schrodinger equation was also 
developed The obtained quantum rate equation can 
properly describe both the many-body Coulomb interac- 
tion and the quantum coherence inside the nano-system. 
However, the rate equation is valid only under large bias 
voltage. 

For mesoscopic transport setup, if we regard the leads 
coupled to the central system as an environment, sim- 
ilar to the consideration for a quantum open system, 
we can perform a perturbative expansion for the tunnel- 
coupling Hamiltonian, and obtain a master equation for 
the reduced density matrix (operator) after averaging the 
states of leads. Indeed, along this line, the perturba- 
tive master equation approach to quantum transport has 
been constructed [p|-p!2[. In most cases, such as the vari- 
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ous quantum open systems as in quantum optics, we find 
that the perturbative master equation up to the 2nd- 
order expansion is good enough. In quantum transport, 
the 2nd-order expansion corresponds to, particularly, a 
sequential tunneling regime. Also, the 2nd-order per- 
turbative master equation covers all the elements of the 
rate equation approach mentioned above. Moreover, the 
number-resolved version of the perturbative master equa- 
tion |]lT|-p!3|, has been demonstrated very successful for 
studying quantum noise and counting statistics in meso- 
scopic transport, including also the large- derivation anal- 
ysis as proposed very recently Q . 

Unfortunately, similar to the rate equation, the 2nd- 
order master equation does not include the tunneling- 
induced level-broadening effect, which implies then a va- 
lidity condition only under large bias voltage. For trans- 
port through Coulomb-island system, on the other hand, 
while the 2nd-order master equation captures well the 
Coulomb staircase phenomenon, it cannot however de- 
scribe the cotunneling process and the nonequilibrium 
Kondo effect. For the latter purposes, higher-order ex- 
pansions of the tunneling Hamiltonian are needed, as pre- 
sented in a variety of scenarios in literature |p|-p!Q|, |l^, ^ . 

In the simpler case of transport through noninteracting 
systems, exact master-equation approach is possible 
In this approach, the non-perturbative Feynman- Vernon 
influence functional theory is employed to derive an exact 
master equation for the reduced density matrix, which is 
valid for general spectral density and for arbitrary tem- 
peratures and bias voltages. Using this exact master 
equation, non-Markovian transport behaviors can be in- 
vestigated. For noninteracting transport under arbitrary 
bias voltage, other exact approaches, being equivalent to 
each other in essence, are also possible. For instance, the 
single electron wavefunction (SEWF) method [^, start- 
ing with a full Schrodinger equation, can evaluate ex- 
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actly the current contribution of an individual electron. 
The total current, in noninter acting case, is thus an in- 
tegration of the individual ones between the (arbitrary) 
voltage window. 



Now we raise a question: is it possible to construct a 
simple master equation, being applicable to interacting 
system and to arbitrary voltages (not suffering from the 
restriction of large bias limit)? In this paper, we make 
an effort aiming to this goal. We notice that, actually, 
the 2nd-order master equation is obtained from the well- 
know Born approximation in the perturbative expansion 
of the tunneling Hamiltonian. The resultant dissipation 
term, in analogy to the quantum dissipative system, cor- 
responds to a self- energy process of tunneling. On the 
other hand, it is well known that in Green's function 
theory, an efficient correction to the self-energy diagram 
under Born approximation is by a self-consistent self- 
energy diagram, under the so-called self- consistent Born 
approximation (SCBA). Therefore, based on this insight, 
for the quantum transport we replace the free (system 
only) Green's function in the 2nd-order self-energy dia- 
gram, with an effective propagator defined by the 2nd- 
order master equation, which favorably includes the tun- 
neling self-energy in it. We found that the effect of this 
simple improvement is remarkable: it can recover not 
only the exact result of noninteracting transport under 
arbitrary voltages, but also the cotunneling and nonequi- 
librium Kondo effect in Coulomb interacting systems. 



The paper is organized as follows. In Sec. II we present 
the main formulation of the master equation approach 
under SCBA, without specifying the concrete systems. 
The presentation in this central part constitutes a num- 
ber of subsections: we first outline in 11(A) the master 
equation approach under Born approximation to quan- 
tum transport, then discuss in 11(B) the basic idea of the 
SCBA which is further implemented in 11(C) to construct 
the improved master equation approach; subsequently, in 
11(D) and 11(E), we consider the steady state and prove 
an exact equivalence with the nGF approach for nonin- 
teracting systems. In Sec. Ill, we continue a challenging 
test for the proposed approach, by applying it to the 
transport through an interacting quantum dot, showing 
a recovery of the nonequilibrium Kondo effect. In Sec. IV 
we apply the approach to an Aharonov-Bohm double-dot 
interferometer, first recovering the exact result of nonin- 
teracting dots, then extending the study of interacting 
double-dots from large bias limit to arbitrary voltage. 
Finally, in Sec. V we summarize the work. As comple- 
mentary materials, we arrange three Appendices (A, B 
and C) for some particular details which are omitted in 
the main text. 



II. FORMULATION UNDER SCBA 

A. Scheme under Born Approximation 

Let us start with a transport setup described by 

H = Hs{ala^) + HB+H'. (1) 

In this Hamiltonian, Hs is for the central system em- 
bedded between two leads which are regarded as a gen- 
eralized bath (or reservoir) and are modeled by Hb = 

T.a=L,RT^fik^c.^^kbl^f.ba^k' And, the coupling between 
the system and the leads is described by the tunnel- 
ing Hamiltonian, = X]c.=L,i? + H.c.). 
Here, aj^ (a^) is the electron creation (annihilation) op- 
erator of the system state "/i", while b^^^j^ (^a/xfe) are the 
creation (annihilation) operators of the leads, associated 
with the labeled specific states. 

Following the conventional treatment of quantum open 
systems, we introduce a collective reservoir (lead) opera- 
tor coupled to the system state "/i", /^^ = t^^kba^k^ 
where a = R. Using it, we can rewrite the tunnel- 
ing Hamiltonian as = Ea=L,i? E/, + H.c.) . 
Then, treating H' perturbatively up to the second-order, 
i.e., under the Born approximation, we formally obtain a 
master equation as |l9|: 

m = -iCp{t) - f dr^^^Ht - t)p{t). (2) 

Here, the reduced density matrix of the system^ p{t), is 
defined by tracing out the reservoir states from the entire 
system-plus-reservoir density matrix prif)-, i-e., pit) = 
TrB[pT(t)]. Two superoperators, the system Liouvillian 
C and the 2nd-order self-energy superoperator Yj^'^\ are 
involved in Eq. (^, with explicit definitions as C{- ■ ■) = 
[HsA-")]. and I](2)(t-r) = (t)a(t, r)^ (r)). Here, 
in E*^^^(t — r) another two superoperators are introduced: 
the Liouvillian £'(•••) = [i^^, (•••)], and the free (sys- 
tem) propagator Q{t^r) = g-^^C^-'r). 

Reexpressing the Liouvillian self-energy superoperator 
in a form in Hilbert-space, the integrand in Eq. (||) has 
four terms, see Appendix A and Fig. 1 where a diagram- 
matic technique on the real-time Keldysh contour is em- 
ployed for an illustrative purpose. Or, in an algebraic 
form, we can obtain an explicit expression for the master 
equation as 

m = -icpit) - J2 { [<,4pH*)] + H-c.}. (3) 

Here, for the sake of brevity, we make a couple of con- 
ventions: a = + and — , a = —a; = aj^ and 

a~ = a^. Also, we denoted A^p^(t) = Y.a^L,R^^^i^p^^)^ 
while A^aiip(t) reads 

4;VW = E f drCl^l(t-r){Qit,T)Kp{T)]}, (4) 
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where the reservoir (lead) correlation functions are de- 
fined simply as C^Ju{t - r) = (/a^^ (t)/^^^ (r))B. As in- 
troduced above, here we make similar convention that 

fLVi'^) = /i/i(^) and f^^\t) = fa^l{t)' The time depen- 
dence of these operators originates from using an inter- 
action picture with respect to the Hamiltonian of leads, 
while the average {• • •)b is over the bath (leads) states. 

It should be noted that, in deriving the above results, 
only the Born approximation was used, but not involv- 
ing the Markovian approximation. The non-Markovian 
feature is manifested by the time non-local form in the 
self-energy terms, Eq. (Q). 

Finally, following Ref. associated with the above 
master equation approach, the transport current simply 
reads 

^"(^) = I {Tr[4+),(t)a^ - A(-),(t)at ] } . (5) 

/i 

Here, the trace is over the states of the central system. 



B. Basic Considerations 

In this subsection, we present the basic considerations 
that can lead us to an improved scheme of master equa- 
tion. Strictly speaking, the 2nd-order master equation 
applies only to transport under large bias voltage. That 
is, the Fermi levels of the electrodes should be consider- 
ably further, at least several times of the level broadening 
r, away from the transport levels of the central system. 
This can be understood by the following simple example, 
say, resonant transport through a single-level quantum 
dot. 

For this simple case, one can easily check that the 2nd- 
order master equation would predict a vanishing cur- 
rent at zero temperature, if the dot level Eq is a lit- 
tle bit higher than jiL {fJ^L > IJ^r)- However, it is well 
known that a full quantum mechanical treatment will 
give a nonzero resonant tunneling current in this situa- 
tion. Similar difficulty appears as well if the dot level 
^0 is narrowly in between the voltage window (i.e., (jll > 
Eq > /hr): no matter how small the Fermi levels are away 
from a full resonant current, / = eTL^R/{TL + F^^), 
will be predicted by the theory. Obviously, these unrea- 
sonable results are caused owing to the neglect of the 
level broadening effect. 

We notice that the tunneling self-energy operator in 
Eq. d), E(2)(^ -r) = {C {t)g {t , r) C (r)) , contains a free 
(system only) Green's function Q{t^r)^ which simply cor- 
responds to, in the case of single- level dot, e-*^o(t-r) 
Then, if we attach e~^^^~^^ to this unitary propagator, 
the level (^o) broadening effect is restored. Insightfully, 
the inserted factor, e"^!^"'^!, represents nothing but a 
self- consistent self-energy correction. That is, it improves 
the self-energy diagram, using a Feynman diagrammatic 
language, from the Born to a self- consistent Born ap- 
proximation. 



In Green's function theory, the correction of the self- 
energy diagram under self-consistent Born approxima- 
tion (SCBA) is typically an efficient technique. In Fig. 
1(A), taking the electron-phonon (or photon) interaction 
as an example, we show the basic idea of such correc- 
tion. There, the single (solid) line is for the free Green's 
function Go of the electron, the double-lines for its full 
Green's function G, and the dashed line for the Green's 
function of the phonon (or photon). The basic idea of 
the self-consistent Born approximation is to replace Go 
in the self-energy diagram by an improved one, as 
schematically shown in Fig. 1(A). It is well known that 
this correction corresponds to a systematic re-summation 
of the perturbation and can greatly improve the results 
in most cases. 

For the 2nd-order master equation, Eq. (^, the 
Keldysh diagrammatic representation of the tunneling 
self-energy term is shown in Fig. 1(B). Analogously, 
within the same spirit of the self-consistent Born approx- 
imation, an improved treatment for the tunneling per- 
turbation is proposed in this work, by replacing the free 
(system only) propagator Q{t^r) with an effective one to 
include the tunneling self-energy in it, as schematically 
shown in Fig. 1(C). The resultant scheme, in the remain- 
ing parts of this work, is to be termed as SCBA mas- 
ter equation (SCBA-ME). In next subsection, we shall 
present the detailed formulation simply based on this 
type of observation. 

C. Improved Scheme under SCBA 

Based on Eq. (||) we formally introduce an evolution 
operator U{t,r), which propagates the state in terms of 
p{t) = U{t^r)p{r). Then, the insight gained above leads 
us to replacing Q{t^ r) with U{t^ r), in the self-energy op- 
erator or more precisely in A^J^p [c.f. Eq. (^]: 

At{t) = Yl fdrCi%{t-r){U{t,r)Kp{r)]}. (6) 

Incorporating this improved quantity into the master 
equation and the transport current, we have 

m = -iCpit) - ^ { [a;,A(;J{t)] + H.C.}, (7) 

/icr 

and 

^"(^) = f E K+^p W«M - Ai-liH] } . (8) 

/i 

We notice that, desirably, Eqs. (0) and (js]) have the 
same compact structures as their counterparts under the 
simple 2nd-order Born approximation. The only dif- 
ference is the replacement of A^p{t) by A^p{t). The 
most obvious consequence of this replacement is that the 
broadening effect^ induced by the tunnehng processes, is 
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FIG. 1: (color online) (A): The self-consistent Born approx- 
imation in Green's function theory, where the free Green's 
function Go is replaced by an effective one, G^^\ The 
dashed curve represents, for instance, the Green's function 
of phonon/photon for an electron-phonon/photon interacting 
system. (B): The 2nd-order tunnehng self-energy diagram, 
Tj^'^\t — t), in the real-time Keldysh contour representation. 
The dashed line represents the reservoir electron Green's func- 
tion. (C): The improved tunneling self-energy diagram under 
the self-consistent Born approximation, in which the free (sys- 
tem only) Green's function Q(t,r) was replaced by the 2nd- 
order effective propagator U{t, r). 



included in the system state evolution in the self-energy 
terms of the master equation. But, not only limited to 
this, it has more implications. For instance, a careful in- 
spection of Fig. 1(C) reveals that this replacement, sig- 
nificantly, accounts for the interplay of the multiple tun- 
neling processes and the Coulomb interactions inside the 
central system. It is well known that such type of inter- 
play is the basic reason for Kondo effect, including partic- 
ularly the nonequilibrium Kondo effect in quantum trans- 
port though Anderson-type impurities. In Sec. Ill we will 
detail an example for this understanding. Moreover, it 
can be proved that cotunneling processes are most natu- 
rally contained in the proposed SCBA-ME scheme. Un- 
der the bias condition of Coulomb blockade, the SCBA- 
ME can recover the cotunneling results given by other 
approaches 

Below we outline a protocol of solving Eq. (0), in fre- 
quency domain, by the use of Laplace transformation: 
p{uj) = L[p{t)] = dte^^^ p{t). Performing the Laplace 
transformation on Eq. (^, we obtain 



i(jjp{(jj) — p(0) = —iCp{(jj) 



K-^)]}. (9) 



where A^i^){lo) = Z^q, with A, 
expressed as 



'iilpiuj) explicitly 



E 



duj' 



r(±)(c.O^(c^±c.O[a^p(a;)]. 



In deriving this result, we used 
L[e^'^'^U{t)] = U{uj ± uj'), and 
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-oo 27r^ ^ ' 
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introduce T)^^u{oj) 
5{uj-ek), f^\uj) = fa{^) and /, 



(10) 

simple relation 
Fourier expan- 

In the latter 



where Ta,^y{uj) 



2^ Xl/c '^cx^ikt*avk 

1 — fa{^) are, respectively, the spectral density function 
of leads and the occupied and unoccupied Fermi func- 
tions. 

Eqs. (^) and (1^) constitute a closed form of coupled 
equations, which allows for straightforward evaluation 
for the solution in frequency domain. However, in doing 
this, we must explicitly identify the propagator U{uj) in 
Eq. ([iq| ) . More specifically, we need to consider the evo- 
lution of pj{t) = U{t^ 't)WZp{^)]' Here, for brevity, we use 
a single index "j" to denote the double ones, "{z/, a}". 
In this context, a special care need to be emphasized. 
Originally, the 2nd-order reduced propagator U was in- 
troduced (in time domain) from the usual propagation 
of a physical state, i.e., p{t) = U{t^to)p{to)' However, 
in Eq. (||) or (p!o|), the quantity being propagated is a^p, 
which is not a physical state. Mathematically, we may 
expect from the well-known experiences that the propa- 
gator must be independent of the initial condition, which 
is, in the present context, the object (initial state) to 
be propagated. In most cases, this statement is true. 
However, the analysis in Appendix A shows that this 
"general" rule, quite unexpectedly, breaks down in our 
case. In Appendix A, we prove that, the evolution rule 
(equation) of pj{t) differs from the one for p{t). To our 
knowledge, this is the first time to recognize this sub- 
tle issue in constructing a master equation. We will see 
that, in the following illustrative examples, this clarifi- 
cation is extremely important - otherwise we cannot get 
the correct results. 

Now we summarize the result derived in Appendix A 
as follows: 

^,{t) = -iCp,{t) - ^ [{a„Al^]} + {alAl^]} 



(11) 



The operators A^^^, in this equation are actually the same 

as A^^^ [c.f. Eq. (^], only replacing p by pj. Out of our 
expectation, we notice a significant difference between 
Eq. (pj]) and Eq. (|): the commutators in the standard 
master equation (|) become now the anti- commutators 
in Eq. (|ll|)! 

The solution of Eq. ([ll| ) determines the propagator 
U{uj) in Eq. (p!Q|). Therefore, we Laplace transform 
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Eq. (y): 

—iujpj{uo) — Pj{0) = —iCpj{uj) — T,{(jo)pj{uj). (12) 

Here, in frequency domain, the self-energy superoperator 
reads 



■ (13) 



In this equation, the shorthand notations a^O = a^O 
and a^O = Oa^, considerably simplify the result. And, 
C^ilv{^^\ the Laplace transformation of C^Ju{t), are re- 
lated with r[^u{uo) through the well known dispersive 
relation: 



-c 



dio' 



2n uj±uj' + iO+ 



ri^l(^)- (14) 



Eqs. (^), (|^) and ( pj]) constitute the main formula- 
tion of the SCBA master equation. In Appendix B we 
prove a relation of the SCBA-ME with an alternative 
approach developed recently, say, the hierachical master 
equations (HMEs) scheme The HMEs, in terms of 
a hierachical equation-of-motion (HEOM) for the den- 
sity operator and a series of auxiliary operators, is de- 
rived by a series of derivatives on the Feynman- Vernon 
influence functional. Since the HMEs scheme largely de- 
pends on a spectral decomposition technique, this type 
of construction makes the numerical implementation very 
tricky and time-consuming. In Appendix B we show that 
the SCBA master equation contains, significantly, the 
dominant contribution (terms) of the HMEs. This iden- 
tification indicates that the SCBA-ME scheme should be 
a promising tool for quantum transports. In addition to 
the advantage with transparent physics, the SCBA-ME 
approach also allows for devising more efficient computa- 
tional schemes, by virtue of its compact and closed- form 
structure. 

Based on Eqs. (|^), (|l^) and ( p!^ ) one can actually 
carry out the time-dependent solution. In practice, from 
Eq. ( p!^ ) , we first calculate and store the reduced evolu- 
tion operator, UiuOrn)^ for a series of frequencies (cj^n)- 
Then, we solve the coupled Eqs. (^ and ( p!o) w ith the 
help oiU{(jJm) to evaluate the integral in Eq. (10|). Using 
the solution (in frequency domain), an inverse Laplace 
transformation then gives the system state p{t) and the 
transient current /(t), simultaneously. In this context, 
we would like to mention that, instead of the Fourier 
expansion for C^v(t)^ we can employ the Fade decom- 
position technique as exploited in Ref. [|o| to improve the 
algorithm efficiency. Since this technique can replace the 
integral duo' in Eq. (^0|) by a discrete sum, favorably, 
it reduces the numbers of UiuOm)' 



D. Steady State 

Based on Eqs. (^), (^ and (|l2|) the steady-state result 
can be more easily obtained than the above transient so- 
lution. For this purpose, let us first consider the integral 
in Eq. (^). Since physically, the time correlation func- 
tion in the integrand of Eq. (H) is nonzero only on finite 
timescale (in Markovian limit which is in fact very short), 
we can then replace p{t) in the integrand by the steady- 
state density matrix p, in the long time limit, t ^ oo. 
After this replacement, we first make a Fourier expansion 
for C^Ju{t — r), then Laplace transform U{t — r)[a^p], 
yielding 



afip 



it 



/OO 7 
^^Ti%i{u;M±u;)[atp]. (15) 



Together with Eq. (12), substituting this result into 
Eq. (|^) , we can then straightforwardly solve for p and 
calculate the steady-state current. 

Now, in order to provide further insight for the pro- 
posed scheme, in the steady-state case, we recast the 
current formula Eq. (H) into a more conventional form. 
To this end, we introduce: (fi/^ui^) = tr [a^piiy(cj)] , 
and (p2pu{^) = tr[a^p2zyM], where pii.(cj) and p2u{o^) 
are the solution of Eq. ( [l2|) , with an initial condition 
Pii/(0) = pal and p2u{^) = cttp^ respectively. Moreover, 
we introduce a matrix notation using, for instance, cpj to 
denote the matrix with elements (fjfj^i^. Then, the steady- 
state current can be expressed in a compact form as 



la = ^Ke 
h 



2^ 



Tr{r,(c^) [Uuj)^{uj)-^^{uj)]], 



(16) 



where ^{uo) = ^i{(jo) + <^2(^)- 

Further simplification is possible, if = AF^^, where 
A is a constant. In this case, Eq. (|l^) can be recast to the 
Landauer-Biittiker type of current formula. That is, the 
current is an integration of tunneling coefficient over the 
bias window: I = ^fRe ^ [h{uj) - fniw)] T{uj). In 
our case, the effective tunneling coefficient reads T{uj) = 
Tr{rLri^(rL + TR)-^Re^{uj)]}. Compared to the re- 
sult of nGF formulation [||, we find that plays a role 
of the retarded Green's function, i.e., cp{uj) = iG^{uj). 
The point is that, the current formula in terms of nGF is 
only a formal expression: it does not hint anything for the 
methods to obtain the various Green's functions. Our cp, 
however, is based on a concrete computational scheme in 
terms of a master equation approach. It is therefore of in- 
terest to investigate to what extent the proposed master 
equation scheme under SCBA is precise. In the follow- 
ing, we will demonstrate: (i)itis exact for noninteracting 
systems, and (ii) it is good enough for interacting systems 
- it can predict, for instance, the nonequilibrium Kondo 
effect. 
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E. Noninteracting System: Recovery of the Exact 
Result under Arbitrary Bias 

Consider, in this subsection, the transport through a 
noninteracting multi- level system, Hs = ^^i,h^j,a^^aj,. 
Straightforwardly, based on Eq. , we obtain the equa- 
tion of motion for (fi as follows: 

-iuj(fi{uj) - (^i(O) = -ih(pi{uj) - iT,o{uj)(pi{uj). (17) 

(pi{0) stand for the initial conditions, (pi^u{0) = 
Tr [a^paj,] and (P2fiu{0) = Tr[a^aJ,p]. The 
tunnel-coupling self-energy XIq reads I]o/ii/(^) = 



r 



2tt ui-Lo' + iO+ ' 



(18) 



Then, based on Eq. ([1 71), summing up (pi{oj) and (p2{oj) 
yields 

cp{uj) =i[uj-h- ^oH] = iCiuj). (19) 

In deriving this result, the cyclic property under trace 
and the anti-commutator, {a^,aj,} = (5^^,, have been 
used. 

Eq. ( p!9| ) is the exact Green's function of transport 
through a noninteracting system. We then conclude that 
our improved master equation approach under SCBA is 
exact for noninteracting transports. Inserting Eq. ([l9|) 
into the steady-state current formula given above, we can 
evaluate the current for arbitrary bias voltage. There- 
fore, very importantly, this achievement is beyond the 
ability of the usual 2nd-order master equation approach, 
which applies only to large bias limit, even for the case 
of noninteracting transport. 



III. TRANSPORT THROUGH AN 
INTERACTING QUANTUM DOT 

In this section, we perform a challenging test for the 
master equation under SCBA. That is, we consider the 
transport through a strongly interacting quantum dot, 
modeled by the well-known Anderson impurity Hamilto- 



0.5 



CD 
O 

c 



0.4 



■D 
C 

o 
O 



0.3 



0.2 



■B=0 



I \ I 



B=0.85 Tesla 

B=1.7 Tesla 

B=2.55 Tesla 

^\ . 
\ i\ 



-0.4 



-0.2 



0.0 

V(mV) 



0.2 



0.4 



FIG. 2: (color online) Kondo peaks in the differential conduc- 
tance, where the magnetic field is introduced to generate Zee- 
man splitting, — en^QUnB. The adopted parameters fol- 
low the experiment in Ref ||2l|]: the temperature ksT = 0.005 
meV, the dot level in the absence of magnetic field eo = —5.2 
meV (here we take the zero-bias Fermi level as energy refer- 
ence), and the on-site Coulomb interaction U = 120 meV. We 
also consider an Lorentzian spectral density for the (identical) 
leads as explained in Appendix C, and assume F = 3.4 meV 
and W = 100 meV. 



Un^ni, = aj^a^ is the number operator and U repre- 
sents the interacting strength. 

Applying the SCBA-ME approach to this system, we 
first note that the correlation function, ci^]/, is diag- 
onal with respect to the spin states, i.e., C^t{t) = 

S^^uCi^\t), and r[^Ju = T^^JSi^u. Then, we specify the 
Hilbert space as follows. Simply, there are four states 
involved in the transport: |0), | t)^ I i) and corre- 
sponding to, respectively, the empty, spin- up, spin-down 
and double occupancy states. Using this state basis, we 
can reexpress the electron operator in terms of a pro- 
jection operator, aj^ = |/i)(0| + (— l)^|(i)(/i|, where the 
conventions ( — 1)^ = 1 and (—1)^ = —1 are implied. 

For a solution of steady state, as Eq. (p!5|), we have 



H. = 



u 

r 



(20) 



A 



(±) 



OO 7 

^T(^J{uj)U{±oj)[ap]. (21) 



The index /i labels the spin up ( ) and spin down ( ) 
states, and Jl stands for the opposite spin orientation, 
denotes the spin-dependent energy level, which may ac- 
count for the Zeeman splitting in the presence of mag- 
netic field (5), e^^i = eo ± QI^lbB. Here eo is the degen- 
erate dot level in the absence of magnetic field; g and 
jiB are, respectively, the Lande-_g factor and the Bohr's 
magneton. In the interaction part, say, the Hubbard term 



Straightforwardly, after some algebra, Z//(±a;)[a^p] can 
be carried out as 



U{u:)[alp] 
U{-u))[a^,p\ 



[A+H|M)(o| + <H(-irM)(M|], 
[A;H|o)(Ml + sH(-iriM){d|], (22) 
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FIG. 3: (color online) Coulomb staircase in the current- 
voltage curve. Inset: the corresponding differential conduc- 



tance (unit: e /h). Here, the results based on Eg. (|24|) are 
plotted against a mean- field Hatree-Fock solution Eg. (p5), 
for the purpose of comparison. Their difference gradually 
vanishes with the increase of temperature, as shown from (a) 
to (d). In the calculation, we consider F^ = F^ = F/2 and 
IJiL — —/^R = eV/2. Taking the zero-bias Fermi level as en- 
ergy reference, we set eo = 2 meV, while assuming F = 0.3 
meV and [7 = 3 meV. 



where 



' n;riHnr^i(c.) ' 



Then, we find the solution of cp{u) as 

z[nr^^(c^)-E^+\c^)](l-n^) 



^(1 - %) 

^ ^ 

(24) 

In this result, = p^^ + pdd and 1 - = + 



Poo- Eq. (|24|) , precisely, coincides with the result by an 
equation-of-motion technique of nonequilibrium Green's 
function 0. 

The solution Eq. (^4|) contains the remarkable nonequi- 
librium Kondo effect, by the terms UTi^{- - -)~^ in the 
denominators of ^{oo). In Fig. 2 we demonstrate numeri- 
cally the Kondo peaks in the differential conductance. At 
high temperatures, the terms UTj^[- • • )~^ would vanish, 
reducing Eq. (p4|) to 



(jj ■ 



^O/i 



U 



^O/i 



(25) 



Here we use ^hf{^)^ instead of <^(cj), to imply the result 
being at the level of a mean-field Hatree-Fock approxi- 
mation. As a matter of fact, Eq. ( p5| ) can be derived from 
the equation-of-motion technique of nGF at lower-order 
cutoff, by using a mean- field approximation How- 
ever, even this simple result goes beyond the scope of the 
2nd-order master equation, by noting that it contains the 
broadening effect. In Fig. 3 we plot the current- voltage re- 
lation based on Eq. ( ^4|) against that from Eq. (^5|). Here 
we emphasize that, in the Coulomb-blockade regime, i.e., 
in the plateau stages, the contribution of cotunneling pro- 
cesses has been automatically included in Eq. (pj). 



IV. TRANSPORT THROUGH DOUBLE DOTS 
IN PARALLEL 



Here, we introduced: Il~^{uj) = cj — e^ — Eo/i(cc;) — E^(c(j), 
and n^^(cj) = (jO — — U — So^(c(j) — Ti^{uj). The self- 
energy Eo/i(a;) is given by Eq. (plsj), while ^^{cjo) is defined 
through 



■I 
I 



dco' 



(^') 



27r uj ■ 



27r UJ — Ed 



■ UJ' 



20+ 



(23) 



In this section we consider a more complex setup which 
consists of double (quantum) dots (DD) connected in 
parallel to the leads. Moreover, the interference loop 
is pierced by a magnetic flux (<l>), similar to the well 
known Aharonov-Bohm interferometer. It was shown by 
the recent studies that this system holds some inter- 
esting properties such as non-analytic current switching 
behavior, anomalous phase shift, and giant fluctuations 
of current. And, all of these are induced by an interplay 
of the inter-dot Coulomb correlation and quantum inter- 
ference. In the present work we extend the study on this 
system from the large bias limit of Ref. |22 , to arbitrary 
bias voltage. 
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Following Ref. |22 , we assume that each dot has only 
one level, £^1(2), involved in the transport. For simplic- 
ity, we also neglect the spin degrees of freedom, whose 
effect, under strong Coulomb interaction, can be easily 
restored by doubling the tunneling rates of each dot with 
the left lead Under these considerations, the DD 

Hamiltonian reads 



Hdd 



E 

/i=l,2 



■ U a|aia2<22 • 



(26) 



Here a\ 2 are the creation operators of the DD states, 
and the second term describes the interdot repulsive 
interaction. As in Eq. (1), we describe the leads by 
= T^ki^Lkb^Lk^Lk + ^Rk^'^Rk^Rk], and their coupling 
to the DD by Ht = Y.^,k {^L^io)^hLk^tR^h^^^a^ +H.c.. 
Usually, the coupling amplitudes t^i^R^ij^ are treated being 
of energy independence. However, in some circumstances 
certain difficulty such as divergence will arise from this 
assumption. To remedy this drawback, one may adopt 
an energy-dependent coupling amplitudes, as to be dis- 
cussed in the following, in the interacting case for more 
details. 

In the absence of magnetic field one can always choose 
the gauge in such a way that all couplings are real. In 
the presence of a magnetic flux <l>, however, the tunneling 
amplitudes between the dots and the leads are in general 
complex. We write tL^R^^ = fL(i?)/ie''^^(^)^, where fL(i?)/i 
is the coupling without the magnetic field. The phases 
are constrained to satisfy + (j^m — 4>L2 — 4>R2 = 
where (j) = 27r<l>/^o- In this work, we only consider the 
case with strong inter-dot Coulomb blockade (^1,2 ^ 
jiiL){two dots cannot be both occupied simultaneously). 



A. Noninteracting Case 

We would like to remind that, the usual 2nd-order 
master equation under Born approximation, is applicable 
only for large bias voltage. To demonstrate the applica- 
bility of our master equation under SCBA for arbitrary 
voltages, we first consider a noninteracting case, = 0, 
which may correspond to the double dots being more 
separated away. 

For this noninteracting system, in principle, a stan- 
dard EOM technique in the framework of nonequilib- 
rium Green's function can carry out an exact solution 
Alternatively, in Ref. [^, using a full Schrodinger 
equation based single-electron-wave-function (SEWF) 
method, simple analytic result was obtained in a rather 
transparent way. The basic idea is that, in the absence of 
Coulomb correlation, the current is contributed indepen- 
dently by each individual electron from the source lead, 
with energy CLk' in the window from fiR to fiL (at zero 
temperature). That is, I{t) = Jl^^ deLk'9Llk'{t)^ where 
is the density-of-states of the left lead. Moreover, 
the SEWF approach allows one to calculate Ik' (t) from a 



time dependent Schrodinger equation for the entire sys- 
tem (dots plus leads). Taking = tLiR)2 = hiR) 
and assuming they are energy independent for simplicity, 
we obtain for a single electron stationary current {t ^ 00) 
as follows [||: 

h' = 4.rn\tL\' [ 1/1 P + 1/2 1' + 2Re(e-'*/i/2)] /D^ 

(27) 

where TLin) = 2ngL{R)\tL(R)\'^ and 



/l,2 = (cLfc' 

D = {2eLk' 



El 



Here s = Ei - E2, T = + rR and r<^ = Tl + Trc''^. 

In large-bias limit, \fiL,R — -E'i,2| ^ T, the bounds of 
the integral / deik' can be extended to infinity, which 
enables to obtain an analytic result for the total current 




TLTRsin' 



^TLTRsin^ 



(28) 



where /q = 2r lTr/T is the resonant current in the ab- 
sence of the magnetic filed. The (/)-dependence in Eq. (|2^ ) 
is an example of the Aharonov-Bohm effect. In Fig. 4 we 
illustrate this Aharonov-Bohm oscillations in large-bias 
limit, and the results for finite bias voltages by numeri- 
cally integrating Eq. (27). 

Now, we consider using the SCBA-ME approach to 
solve this same problem. To demonstrate the general al- 
gorithm, which obviously has an implication of applying 
to other systems, we would like to perform a direct nu- 
merical calculation, rather than making efforts to extract 
out any analytic results. In the noninteracting case, we 
apply the state basis |0) = |00), |1) = |10), |2) = |01) 
and 1 3) = |11), to convert the density operator and other 
operators into matrix forms. Here 1 0(1)0(1)) corresponds 
to vacant ("0") or occupied ("1") dot state. Then, com- 
bining Eq. ( p!5|) together with its auxiliary equation ( p!^ ) 
and Eq. (Q) in the long time limit {p = 0), we can find 
the steady-state solution, for both p and A^J^lp. Finally, 

inserting the obtained A^^^^p into Eq. (^), we are ready 
to calculate the steady-state current in terms of a simple 
trace manipulation. 

Figure 4 shows a full agreement between the result 
from the SCBA-ME approach and the exact solution 
Eq. (I27I), for arbitrary bias voltages. This achievement 
desirably goes beyond the 2nd-order master equation ap- 
proach, as employed in Re f. [[22| , where the result be- 
comes consistent with Eq. (^7D only in large bias limit. 
Concerning the interference pattern in Fig. 4, we would 
like to remark that in general the Onsager relation would 
lock the current (or conductance) peaks at = 27rn for 
any two-terminal linear response transport |23 , either 



noninteracting or interacting. Yet, the Onsager relation 
cannot apply to finite voltage transport. While here in 
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Fig. 4 the current peaks are still at = 27rn for finite 
voltages, we found that this is not the case for interact- 
ing double dots. There, the phase- locking behavior as 
preserved in Fig. 4 breaks down p2| . 



□ v=r < v=4r o v=20r 




— I 1 1 1 \ 1 1 1 1 — 

-2-1012 

FIG. 4: (color online) Aharonov-Bohm oscillations of current 
through noninteracting double dots. We assume a symmetric 
setup with Tl ^Tr = P. lo = eFLTR/irL + F^) = eF/2 
is the large-bias-limit current in the absence of the magnetic 
flux. Also, we set ijll — —/J^r — Ei — —E2 — ^5 
A = O.IF. We find that excellent agreement between the ex- 
act solution (solid lines) and the numerical SCBA-ME results 
(symbol curves) is achieved under arbitrary bias voltages, as 
exemplified by the three values. 



an unexpected switching behavior holds only for large 
bias voltage. 

For finite small voltages, unfortunately, either the 2nd- 
order master equation or the SEWF approach is no longer 
apphcable. Now, beyond Ref. [|2|, we apply the SCBA- 
ME approach to this unexploited regime. The calculation 
scheme is straightforward, simply along the line of solv- 
ing the noninteracting DD, as described in the previous 
subsection. However, there is an important issue need- 
ing care in this context. First, notice that our SCBA- 
ME is constructed in terms of a time-nonlocal version, 
see Eqs. (^ and (|^), where p(r) enters the integrand of 
Jo dT[- • •]. For this time- nonlocal version of master equa- 
tion, we know that even at the usual 2nd-order level, 
the lower bound of the integral dT[- - -] would re- 
sult in a "principle- value" contribution, or, renormal- 
ization effect. For our interacting DD interferometer, 
as analyzed in detail in Ref. [^, the current switch- 
ing behavior is caused by the occupation of a trans- 
formed dot state. On the other hand, however, we can 
verify that the "principle- value" integral will induce an 
effective coupling between the transformed dot states, 

^12 = r dcLkP^ + dcLk ^7#r7n- Then, if 

JfJ^L ^'^ E2-eLk J-00 ^'^ eLk-{Ei-\-U) ' 

we adopt the usual wide-band-limit model for the leads 
(i.e., constant density-of-states and energy-independent 
tLk)^ we will find that the integral can become even di- 
vergent! 



B. Interacting Case 

When the double dots are put closer, strong Coulomb 
repulsion among the inter dots may lead to a reduced 
Hilbert space expanded only by {|00), |01), |10)}. In this 
large U hmit, as well as under large bias voltage {U 
I Ml, i? — ^1,2! ^ F), an analytic result for the steady-state 
current was carried out in Ref. p2|- 



m = Ic 



£2 + Ic (2rR sin^ I - e sine/)) ' 



(29) 



where Ic = 2TlTji/ {2Tl + F^^) is the current in the ab- 
sence of the magnetic flux. Comparing Eq. ( p9| ) with 
Eq. ( p^ ) for £ 7^ we find that both currents display 
the Aharonov-Bohm oscillations. Nevertheless, Eq. ( |29| ) 
shows an asymmetric behavior with respect to the mag- 
netic flux {(f)). This can be verified by finding, from 
Eq. (I29I), the current peaks not at = 27rn. This 
phase-blocking violation is an interesting consequence of 
Coulomb correlation in the double-path interferometer, 
differing from the noninteracting double dots as discussed 
above. More strikingly, when £ ^ 0, Eq. ( p9|) predicts 
a non-analytic current switching behavior. That is, one 
finds that / = /c for = 27rn, but / = for any other 
value of (j). However, as to be shown in the following, such 



This difficulty can be eliminated by applying the more 
realistic tunneling Hamiltonian with energy-dependent 
amplitudes (tLk)- Unfortunately, to our knowledge, so far 
the specific behavior of tLk with the energy is not well 
clarified. Therefore, based on certain intuition, in this 
work we adopt an ansatz, \tLk\^ = |tope"^''''^"^^^'/'^', 
where AE is the DD energy change associated with an 
electron jump, into the dots or out to the leads. Incor- 
porating this ansatz into Eq. (|5|), we attach to T^u{^) 
with a decay factor, e~(^~^"^") , where the energy dif- 
ference Emn = Em — En corrcspouds to the DD states 
|m) and |n), in association with the matrix element 

ran • 



Figure 5 displays the current of the interacting DD 
interferometer with the magnetic flux, under finite bias 
voltages. We would like to mention that, for arbitrary 
voltages, the current pattern for e = 0, aligned DD levels, 
will keep the current peaks at = 27rn. For £ 7^ 0, how- 
ever, the current peaks will have an e-dependent shift, 
similar as in the large bias limit. We see that, in Fig. 
5, the current under finite voltage has a relatively sharp 
change in the vicinity of ^ = 27rn, which is an omen for 
the current switching in large bias limit. With the in- 
crease of the bias voltage, the switching behavior is to be 
more pronounced. 
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Appendix A: Derivation of Eq. (11) 



FIG. 5: (color online) Correlation and interference induced 
switching behavior of current through the interacting double 
dots. Parameters and conventions are the same as in Fig. 4, 
except that a smaller e — COIF is assumed in order to result 
in the current switching phenomenon in large bias limit. We 
assume the width parameter a — lOF in the coupling- strength 
ansatz. 



V. SUMMARY 

We have proposed an improved master equation ap- 
proach to quantum transport, by implementing a general- 
ization from the Born to self-consistent Born approxima- 
tion. Being in reasonable agreement with other scenar- 
ios of the self-consistent Born approximation, we found 
that the proposed scheme can give satisfactory results. 
For instance, it can recover not only the exact result of 
noninteracting transport under arbitrary voltages, but 
also the cotunneling as well as the challenging nonequilib- 
rium Kondo effect in Coulomb interacting systems. Very 
importantly, this achievement is completely beyond the 
scope of the usual 2nd-order master equation under Born 
approximation. 

As a final remark, we notice that, while the correction 
from the Born to self-consistent Born approximation is 
an efficient and well-known technique in Green's function 
theory, similar correction in master equation is lacking so 
far among the various efforts. Since it is impossible to 
construct an exact and closed form of master equations 
in general, the scheme proposed in this work, containing 
the dominant contribution, is extremely desirable. 



In this appendix, we present the derivation of Eq. (O), 
i.e., the equation-of-motion of pj{t) = Z//(t, r)[a^p(r)]. 
Similar to deriving the usual 2nd-order master equation, 
we start with a 2nd-order expansion for the tunneling 
Hamiltonian H\ and formally obtain the same equation 
as Eq. (|): 

^j{t) = -iCpj{t)- [ dTTTB[C\t)g{t,T)C\r)pT{r)]. 
Jo 

(Al) 

The various superoperators in this equation have the 
same meaning as in Eq. (||). Explicitly, the Liouvillian 
self-energy superoperator can be reexpressed in Hilbert 
space as 

£'{t)g{t,T)C'{T)pT{T) 

= H'{t)G{t,T)H'{T)pT{r)GHt,r) 

-G{t,T)H'{T)pT{T)G\t,T)H'{t) 
- H'{t)G{t,T)pT{T)H'{T)G^{t,T) 
+ G{t,T)pT{T)H'{T)G\t,T)H'{t) 

= 1 -II -III + IV. (A2) 

Applying the Bom approximation, Pt{t) ~ PBPj{T), we 
further obtain 

TtbII] = TVB[H'{t)G{t,T)H'{T)pBPj{T)G\t,T)] 

= T.{^lt\t-T)a^g{t,T)[alpj{T)] 

+ C\,-\t - T)alg{t,T) [a,Pj{T)\ }, (A3a) 

TiBill] = TrB[G{t,T)H'{T)pBPj{T)GHt,T)H'{t)] 
= -T.{^ltKt- T)g{t,T) [atp,{T)]a, 

+ Cl-Jit - T)g{t,T) [a,p,{T)]al}, (A3b) 

TrBlIII] = TVB[H'{t)G{t,T)pBPj{T)H'{T)GHt,T)] 

= - E {^^^^*(^ - [pj{T)al] 



■ClV*{t-T)alg{t,T)[pj{T)a^]], 
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Ttb[IV] = TTB[G{t,r)pBPj{r)H\r)GHt,r)H\t)] 

IJ,U 

+ Cl+>{t-T)g{t,T)[pj{T)a,]al}. 

(A3d) 
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Substituting these results into Eq. (|Al|), a more compact 
notation leads to Eq. ([ll|). 

It is of crucial importance to note that, in the above 
Ttb[II] and Tr^ [///], an extra sign change occurs when 
we arrange the bath operators to the two sides of pB- 
We explain this issue in more detail as follows. Con- 
sider, for instance, TTB{[atfu{r)][pBpj{r)][fl{t)a^]}. In 
order to utilize the cyclic property under Tr^[---], i.e., 

TTB[fu{r)pBfl{t)] = clV{t-r), we have to move /t(t), 
crossing Pj{r), to the right side oi pB- Since Pj(r) con- 
tains a Fermi operator, aj = a^, the above movement 
of /^(t) will cause an additional minus sign, according 
to the Fermi-Dirac anticommutative relation between aj 
and Similar consideration and result apply also to 
all the other terms of Ttb[II] and Tr^ [///]. As a con- 
sequence, the resultant extra sign of minus, very impor- 
tantly, alters the commutators in the usual master equa- 
tion, e.g., Eq. (H), to the anti- commutators in Eq. (pTj). 



Appendix B: Relation with the equation-of-motion 
approach 

The proposed SCBA-ME is based on an insight from 
the Feynman's diagrammatic technique. In Green's 
function theory, it is well known that there exists 
some relation between the diagrammatic technique and 
the equation-of-motion approach, despite the lack of 
a general equivalence between them. Rather than 
the equation-of-motion for Green's functions, a hier- 
achical master equations (HMEs) scheme for the den- 
sity operator and a series of auxiliary operators were 
constructed [|l^, based on the Feynman- Vernon influence 
functional theory for quantum open system. In this ap- 
pendix, we extract out an intrinsic relation between the 
present SCBA-ME approach and the HMEs scheme. For 
this purpose, we introduce first the notation used in Ref. 
^5l: 



(Bla) 
(Bib) 



Here w4a^p(cj,t) is defined from ^i^p(t) = 
/ ^A^a^p{uj,t). Then, from Eq. (|) we have 

Aiti^.t) = ^ Tdrri^Uc^je— (^--) {U{t,r)Kp{r)]} , 

(B2) 

Compared to Ref. ||, we find pi^j(t) = / (cj, t), 

which is nothing but the first tier auxiliary density oper- 
ator introduced there. 

Now, consider the equation-of-motion of 



U{t^r)[a^p{r)]. Formally, from Eq. (|l 1\) , we have 

dtU{t,to) ^ -i€ll{t,to) - i I '^{a%U^^^{^,t)}, 

dtUi^J{w,t) = -i{C - auj)Uil\uj,t) - iCi^^4wMt,to). 

(B3) 

Here, we have further introduced 

Ci;l(c^)0 = E K,A^)<0 ± Tl^^iw)da:] , (B4) 

and 

Ui-J{u;,t) = -i rdTe-*(^-'^-)(*--)ci;)+HZi(r,to). 

_ (B5) 
With these identifications, Eq. (Q) can be recast to the 
following set of coupled equations: 

m = -iCp{t) - i E / ^ [<' -^"M . (B6a) 



-iCi^)4u;')<pi-^{u;,t). (B6c) 



Eqs. (B6), favorably, recover the HMEs truncated until 
the second tier (ntier = 2) [|l^, up to an unimportant 
difference in Eq. ( |B6cD . That is, the HMEs with a cutoff 

at ntier = 2 has one more term, C^^_^{uj)(j)^^,^,{uj\t). We 
have numerically examined, however, that this term has 
negligible effect. 

In good agreement with other SCBA scenarios in broad 
physical contexts, we find that the SCBA-ME can give 
satisfactory results, not only for noninteracting trans- 
ports, but also for highly correlated processes. For in- 
stance, it can predict the nonequilibrium Kondo effect 
[||, while this challenging effect is completely beyond 
the scope of the second-order master equation, as well 
as some improved approach We have also checked 
that the proposed SCBA-ME contains all the cotunnel- 
ing processes in Coulomb blockade systems |10 . On the 
other hand, in applying the HMEs, it was found that the 
dominant contributions are usually obtained by a cutoff 
at ntier = 2, while the pure numerical convergence may be 
achieved by including more tiers of the HMEs, at a price 
of much more expensive computations. The above facts, 
combined together, then suggest that the SCBA-ME is a 
reliable quantum transport approach. Therefore, based 
on the SCBA-ME, the above identification of Eqs. (p6|) 
actually helps us extract out the dominant tiers. That 
is, the first two tiers of the HMEs would dominate the 
numerical contribution, as well as the underlying physics. 
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Appendix C: Lorentzian Spectrum 



pressed as 



In most cases for a system coupled to a continuum, 
a Lorentzian type of spectral density function is more 
reasonable than a constant. In quantum transport, for 
coupling with the leads, we therefore assume p5[ 



(CI) 



In certain sense, this form corresponds to a half-occupied 
conduction band for each lead, which fixes the center 
of the Lorentzian function at the Fermi level (chemical 
potential), fia^ of the lead. In Eq. (pT[), Wa characterizes 
the bandwidth of the a-th lead. Quite naturally, the 
constant coupling rate can be recovered from Eq. (CI) by 
assuming Wa oo, yielding Tai^u{(^) = ^ai^w 

Corresponding to the Lorentzian spectral density 
Eq. (CI), the correlation function of Eq. ( [T^ ) can be ex- 



(C2) 



The second quantity is related to the first one through 
the well-known dispersive relation as 



-( 



^^_r(±) H 



Re 



1 /?(c.-Ma) 

2 27r 



'1 ^\ 

2t, r Wo, 



(C3) 



where V denotes the principle value of the integral, and 
^{x) is the digamma function. 
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